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Abstract 

The assumption that a solution to the Einstein equations is static (or stationary) very 
strongly constrains the asymptotic behaviour of the metric. It is shown that one need 
only impose very weak differentiability and decay conditions a priori on the metric for the 
field equations to force the metric to be analytic near infinity and to have the standard 
Schwarzschildian falloff. 
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I: Introduction 

It has long been known that the asymptotically flat static and stationary solutions of 
the Einstein equations are analytic at spacelike infinity (Beig and Simon 1980, 1981). This 
should come as no real surprise, because they are the special solutions without gravitational 
radiation and so therefore the far- field should be dominated by the Newtonian potential(s). 
Further, Beig and Simon (1980, 1981) have shown that all static and stationary solutions 
are Schwarzschildian at infinity. 

In this article, we wish to improve on this result by weakening (in a significant way) 
the assumptions made, yet recovering the Beig and Simon results, both about analyticity 
and about falloff. 

One assumption made by Beig and Simon was that the given (static or stationary) 
metric went flat like 1/r (in some asymptotically cartesian coordinate system). It has been 
shown that 'reasonable' solutions to the Einstein equations exist which fall off much slower 
at infinity (Christodoulou and 6 Murchadha 1981). For example, the condition for finite 
total energy is that the metric go flat faster than r~i (Bartnik 1986, Chrusciel 1986a,b, 
O Murchadha 1986). With such a fall-off, it can be shown that the energy is very well 
behaved, it is positive and is a component of a well-defined timelike energy-momentum 
vector which transforms correctly under Lorentz transformations. Christodoulou and O 
Murchadha (1981) even discussed solutions which went flat essentially arbitrarily slowly. 
Such solutions, despite the fact that their masses were infinite, could reasonably be de- 
scribed as being asymptotically flat with complete spacelike infinities. 

One result of this paper is that all such slow decay at infinity must be due to the 
existence of gravitational radiation near spacelike infinity. In particular, we show that if 
we assume that the solution has slow falloff near infinity and that it is simultaneously 
either static or stationary then the field equations force the usual £ falloff on us and we 
recover the standard Beig and Simon Schwarzschildian result. 

In addition to allowing slow decay at spacelike infinity, we also assume that the metric 
belongs to some Sobolev space and so need only be weakly differentiable, rather than 
belonging to any classical space. Again, the field equations force up the differentiability 
and we can finally conclude that the metric must be analytic. 

In this paper, we are careful to impose conditions only on the metric in the physical 
spacetime. We then use the field equations to show that other (unphysical) quantities 
are suitably smooth as we need them. We take a straightforward 3 + 1 approach to the 
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problem and write the static and stationary field equations as conditions on the Cauchy 
data for the gravitational field (Arnowitt, Deser and Misner 1962 (ADM)). 



II: Asymptotically fiat static and stationary solutions 

We will consistently adopt an initial value approach in this article. This means that 
we assume we are given initial data for the gravitational field on a single spacelike slice. 
The Einstein evolution equations now determine how this data evolves. This evolution 
then can be tested to see whether the spacetime in either static or stationary. We are only 
interested in working near infinity so we will assume that the data is given on a region 
equivalent to R 3 with a ball cut out. 

Initial data for the vacuum gravitational field consists of a pair (5^, 7r lJ ) (ADM), 
where gij is an asymptotically flat Riemannian three-metric on the given manifold and 
7T* J is a symmetric three-tensor density. These are not independent, they must satisfy the 
constraints 

^R= 1 -[^\ ah - 1 {g ah ^f] (2.1) 

9 2 

ViTT^ = (2.2) 

where is the three scalar curvature and V is the covariant derivative compatible with 
gij . The rest of the Einstein equations form the evolution equations 

^ = 2Ng- 1 * [n l3 - l -(g ab n ab )g l3 ] + N Mj + N j;i (2.3) 



^ = -Ngl[B*1 - \jiR] + VW° 6 " \{9a^ ab f] ~ 2Ng^[^ m 

- Iq^V- 7 ] + 9kN' ij - g lJ V 2 N) + (^N m ). m - Nf m * m * - N^ im (2A) 

The extra scalar and vector (AT, N l ) (the lapse and shift) represent the arbitrariness in the 
choice of "time" off the hypersurface. That is 

= iVn M + AT M , n^N^ = 0, (2.5) 

where is the time translation vector and n M is the unit (timelike) normal to the hyper- 
surface; is really a three-vector, it lies entirely in the given slice. 
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Now we wish to consider the case where the given data represents a slice through a 
spacetime which has a timelike Killing vector /c M (at least near infinity). This implies that 
for the given (g, it) there exists a particular choice of (N, N l ) which, when substituted into 
(2.3) and (2.4), gives 

9g - - - 0. (2.6) 



dt dt 

This special choice of lapse and shift are those which correspond to using the killing vector 
as "time" direction. In other words, we split up as suggested by eq.(2.5), to generate 
the desired (N,N*), 

= Nn 11 + iV M . (2.7) 

Let us further assume that the spacetime is static. This means that we can assume that 
the killing vector is hypersurface orthogonal. Let us choose our three-slice to be the 
hypersurface that is perpendicular to the killing vector. This means that the killing vector 
is parallel to n M and the shift iV M will be zero. To repeat: for a static spacetime we must 
have an initial data set (g, it) (satisfying the constraints) and a choice of lapse N (with 
N l = 0) such that (2.6) holds. N is just the length of the killing vector. 
In the light of this (2.3) reads 

N[n lJ -^(g ab n ab )g lJ }=0. (2.8) 

Since we must have N ^ (otherwise the Killing vector vanishes) we immediately get 

n lJ = 0. (2.9) 

Obviously, (2.2) is trivially satisfied. The Hamiltonian constraint (Eq. (2.1)) gives us 

{3) R = 0. (2.10) 

When (2.9) and (2.10) are substituted into (2.4) we get 

= -NR ij + (N' lj - g tj V 2 N). (2.11) 

The trace of (2.11) gives 

= -N( 3) R- 2V 2 iV. (2.12) 
But since ^.R = (Eq. (2.10)) we get 

V 2 iV = 0. (2.13) 



Finally, when (2.13) is substituted back into (2.11) we find that 



NRij = N-ij. (2.14) 

Therefore a static solution to the Einstein field equations is equivalent to having a Rie- 
mannian three- metric g and a scalar function N which satisfy (2.14), together with either 
(2.13) or (2.10). We need only one or the other because the trace of (2.14) gives 

N {3) R = V 2 N. (2.15) 

Having obtained the static field equations in the form (2.10), (2.13), (2.14) we wish 
to manipulate them further. We wish to use the lapse N as a conformal factor to simplify 
the equations. In particular we define a function u by 

e u = N, (2.16) 

and transform the metric by the rule 

g ab = N 2 g ab = e 2u g ab . (2.17) 
A standard calculation now shows 

Rji = Rji - u.ji - g 3 iW 2 u + u.ju-i - (Vm) 2 ^. (2.18) 



It is easy to show 



— U ;l j + U ;l U ; j, 



V 2 iV = e u [V 2 M + (V«) 2 ], 

V 2 u = e- 2u [V 2 u + (Vu) 2 }. (1.19) 
Therefore (2.13), (2.14) corresponds to 

Rij = +2U;iU;j, (2.20) 

and 

V 2 -u = 0. (2.21) 

Eq.(2.20) and Eq.(2.21) express the static field equations in the form that is most conve- 
nient for further analysis. 
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The equations do not simplify so dramatically in the stationary case because the 
Killing vector will not be hypersurface orthogonal. However, some significant simplification 
can still be achieved. We can no longer force the conjugate momentum (n^) to vanish 
entirely but we can always choose our slice so that gijn^ = 0, at least near infinity (Bartnik, 
Chrusciel and 6 Murchadha 1990) (which is all we are interested in). 

Taking the trace of Eq.(2.3) gives 

-Ng-^g ab n ab + 2V a N a = 0, (2.22) 

which immediately gives 

V a N a = 0. (2.23) 

The trace of Eq.(2.4) gives 

lNg$R - ^Ng-K ab n ab - 2g^V 2 N - 2N a:b n ab = 0. (2.24) 
Using the Hamiltonian constraint (2.1) we can eliminate the first two terms to give 

^V 2 AT = -iV a:6 7r a& (2.25) 

If we multiply (2.3) by 7r lJ we get 

Ng-K ab n ab = -N a:b n ab . (2.26) 

If we now substitute (2.26) into (2.25) to give 

V 2 iV = Ng-h ab n ab = NR. (2.27) 

This is the well-known equation for the choice of lapse which preserves the maximal slicing 
(gab^ ab = 0) condition. Finally, we can substitute (2.27) into (2.4) to simplify it and give 

Ng^R ab = g^N, ab - 2Ng-^^n mb + (n ab N m ), m - N a , m ^ - N b , m n^. (2.28) 



III Weakly decaying static solutions 

Following from Section II, we seek a metric g and a function N satisfying 

NR ab = N. ab (3.1) 

V 2 iV= (3) J R = (3.2) 

We now need to define suitable function spaces in which to define them. We choose to use 
weighted Sobolev spaces. Let us define: 
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Definition 1. a = (1 + |x 2 |)i 

Definition 2. H s s( 3 R/ B), s G N, 6 G dt, is the class of all functions f on 3 R/B possess- 
ing weak derivatives up to order s such that a~ 5 ~^ + ^D a f , for each \ct\ < s, is square 
integrable on 3 R/B. 

This forms a Hilbert space with norm 




Note: This definition of a weighted Sobolev space is the standard one (as in Choquet- 
Bruhat and Christodoulou 1981) except that (following Bartnik 1986) we have subtracted 
3/2 from the old definition of 5 and changed sign. Now a function in H s § has (more or 
less) classical blowup like r +s . 

One standard result (Choquet-Bruhat and Christodoulou 1981) we will frequently use 
is the: 

Multiplication Theorem. Pointwise multiplication is a continuous bilinear map 

{hi h) - >• fih 

H Sl ,s 1 x H S2y s 2 — > H SjS 

if si, s 2 > s, s < si + s 2 - 3/2, 5 > Si + S 2 . 

We wish to consider static solutions that are asymptotically flat in the sense that the 
metric (in some coordinate frame) approaches the flat cartesian metric at infinity. We also 
assume that the killing vector approaches the standard time-translation killing vector in 
the same frame. In other words, we wish to consider solutions to (3.1), (3.2), (3.3) on a 
manifold with topology 3 R/B, having the following properties: 

g - 5 G H SjS , s > 3 , 5 < (3.4a) 

N -> 1 at oo . (3.46) 

Condition (3.4a) on the metric says that it approaches the cartesian metric 5 at infinity, 
but we do not care how slowly; the condition (3.4b) on iV (the lapse, the length of the 
Killing vector) says that it goes to a time-translation at infinity. Standard theorems on 
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the laplacian (see e.g. 6 Murchadha (1986)) allow us to deduce immediately from (3.2) 
that N must satisfy 

M 

N - 1 + — e # s+ i, 5 _i , (3.5) 

where M is some constant. 

We know that N — 1 satisfies V 2 (iV — 1) = and when one explicitly writes this out it 
becomes clear that the leading part of N — 1 must be a harmonic function of the flat space 
laplacian. The residue then falls off faster than 1/r. Note also that N (again because it 
satisfies the laplacian) is one degree more differentiable than g. 

Now we know that N is smooth and well-behaved [in particular we need that N is 
bounded away from zero in a neighbourhood of infinity, we can use the standard embedding 
theorem (Choquet-Bruhat and Christodoulou 1981) to prove this], we can use it as a 
conformal factor, g = N 2 g = e 2u g. g can be shown to be a riemannian metric which 
satisfies 

~g - 5 e H s , s , s > 3, 5 < , (3.6) 

outside some ball B' . 

The field equations take the form ((2.20), (2.21)) 

Rij = +2u-iU.j, (3.7) 

and 

V 2 u = 0, u -> at oo . (3.8) 
again, the standard theorem for the laplacian on weighted Sobolev spaces gives 

u - — e H s+lyS -i , (3.9) 

for some constant A. (This also follows directly from (3.5) and u = In AT.) 

When this is substituted into (3.7) we see that the Ricci curvature of the conformally 
transformed space falls off like 1/r 4 . We actually get 

Rij - 2A 2 ^f e H Sj5 -a, (3.10) 

where A is the same constant as in (3.9). Since in three dimensions, Ricci is equivalent 
to Riemann, this means that the Riemann curvature also falls off like 1/r 4 . This, in 
turn means that the space must be flat to order 1/r 2 . When we return to the physical 
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space by multiplying by N 2 we introduce a 1/r term into the metric but this term is a 
pure conformal factor of the form 2M/r. Therefore the physical three-space must be pure 
Schwarzschild to order 1/r, g ab = (1 + 2M/r)S ab + 0(l/r 2 ). 

This argument can be made more precise by introducing harmonic coordinates in the 
conformal space. We seek three functions 0i, 02, 03 satisfying 

V 2 0i = , 0i -> x at oo , (3.11) 

and similarly for 2 and 03. Such functions exist and belong to H s+ i^. We use these as 
our new coordinates near infinity (the details can be seen, for example, in O Murchadha 
1986) and the transformed metric g' will belong to the same space as g, g' — 5 G H Sj s- 
Writing h ab = g' ab - 5 ab , we get 

R ab = V 2 ~h ab + IT. (3.12) 

This is the standard simplification of the Ricci tensor that occurs on using harmonic 
coordinates. V is essentially just the first derivative of the metric so it is easy to show 
r G H s -i y s-i- The multiplication theorem can now be used to show IT G H a -i^5~2- This 
is where we need to use s > 3, otherwise the derivatives go bad. We rewrire (3.12) as 

V 2 h ab = R ab - IT. (3.13) 

Since we know, from (3.10), that the Ricci tensor is both smoother and has faster falloff 
than IT, we see that (3.13) forces h ab into H S+ \ 2S- Therefore we improve both the 
differentiability and the falloff of the metric. We iterate this argument / times where / 
satisfies 2 l 5 < —1. At this point we find that IT (classically) falls off faster than 1/r 3 . At 
this stage when we invert the laplacian (see O Murchadha 1986) we pick up a harmonic 
function of the flat-space laplacian. In other words we have 

h ab - ^ G tf, a ,_ e _i, (3.14) 
r 

where A ab are six (apparently independent) constants. However, we also require that the 
metric be harmonic, i.e., 2h ab , b — h bb ^ a = 0. When we apply these conditions to (3.14) 
we get a set of relationships between the A ab s which cannot be satisfied. Therefore they 
must all vanish and h ab must fall off faster than 1/r. We can now continue the iteration 
and show that h ab must fall off like 1/r 2 . At this stage we have improved the falloff and 
differentiability sufficiently so that we can appeal to Beig and Simon and get that h ab must 
be analytic. 

The content of this Section can therefore be summarized as 
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Theorem. Given an asymptotically Hat Riemannian three-metric g — 5 G H 3j §, S < 0, 
which is the metric on a slice transverse to the killing vector in a static space-time which 
is vacuum near infinity, and given that the length of the killing vector goes to a constant 
at infinity, it follows that the metric must be analytic and Schwarzschildian at infinity. 

Remarks: (i) We need three weak derivatives because we have to estimate terms like 
{9ab,c) 2 and have them smoother than gij^i so as to improve the differentiability. If we were 
willing to consider non-integral differentiability what we require is x — 2 < 2 (x — 1) — 3/2. 
This implies x > 3/2. In particular s = 2 will work. 

(ii) The 5 < condition classically means that we are willing to consider 
any metric which decays to flat space like l/r e for any e. Obviously one would like to 
replace this with 'going flat' and not require any kind of power-law decay. It is difficult to 
see how this might be achieved; none of the present battery of weighted spaces (classical, 
Holder . . . ) seem suitable. 

(iii) The condition N — > 1 at infinity can be significantly relaxed. All we 
need to do is to eliminate the harmonic functions of the laplacian. Therefore all we need 
is that N grows slower than r at infinity. The question of the behaviour of the solution if 
one places no restriction at all on N must be analyzed separately. 



IV: Weakly decaying stationary solutions 

As we showed in Section II, the stationary equations are equivalent to the existence 
of a Riemannian metric g, a tracefree symmetric tensor P (we define P ab = g~^Tx ab to get 
rid of the determinants of g), a scalar N (the lapse, the length of the part of the killing 
vector orthogonal to the slice) and a vector N l (the shift, the part of the killing vector in 
the slice) satisfying 

^R = P ab P ab . (4.1) 
V a P ah = 0. (4.2) 
V 2 A = NR. (4.3) 
V 2 A = -N a . b P ab . (4.4) 
2NP ab + N a , b + N b , a = 0. (4.5) 
g ab P ab = 0. (4.6) 
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V a N a = 0. (4.7) 

NR ab = N. ab - 2NP™P bm + P ab ., m N m - N avm P? - N b , m P™ • (4.8) 

As in Section III we want to make some (weak) assumptions about the asymptotic be- 
haviour of the metric, the extrinsic curvature and of the killing vector. We wish to assume 

9ab - S ab e H 3jS , (4.9) 

N -> 1, N a -> at co. (4.10) 
Equation (4.9) gives us that G This essentially forces (from eqn.(4.1)) 

P ab eH 2A _ v (4.11) 

We assume that this holds. Equation (4.1) also tells us that > and this means that 
N (from eqn.(4.3)) satisfies a maximum principle from which we find that (Theorem 3.5, 
Christodoulou and 6 Murchadha 1981) 

N-leH 3iS . (4.12) 

Taking the divergence of the tracefree part of (4.5) gives 

-2N. a P b a = [N a , b + N b . a - hj c N c g ab r. (4.13) 

The multiplication theorem, together with eqns.(4.11) and (4.12), gives us that 2N. a P b a e 
H 2 M_2- The right- hand-side of (4.13) is the divergence of the conformal killing form of 
the vector iV a . This is a very nice second order linear elliptic operator. We immediately 
get (Lemma 3.1, Christodoulou and O Murchadha 1981) 

N a eH 3 ^_. (4.14) 

We can substitute (4.11) and (4.14) into (4.4) to get V 2 iV e -#2,25-2- In turn this gives 

N - 1 e #4,25. (4.15) 

Comparing with (4.12), we see that we have improved both the smoothness and the falloff 
of N. When (4.15) is substituted into (4.13) we now get that 2N. a P b a e H 2 ss__ 2 . In turn 
this gives 

N a e H 3 5s. (4.16) 
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We can iterate this procedure and keep improving the falloff of both N and N a . There 
exists an integer n such that nS < — 1. We get V 2 iV e H 2 , n 5-2- Then we find 

N-l + C/r eH 4M , (4.17) 

where C is some constant. We also find that N a falls off like 1/r. Thus we can show that 
the killing vector must be quite regular at infinity. 

We now need to show that the metric and extrinsic curvature have similar properties. 
The easiest way to proceed is to imitate Beig and Simon. We have shown that N and N a 
are well-behaved at infinity. Therefore we see that A = N 2 — N a N a is also regular, so that 

X-l + C/r e# 3 ,-i-e, (4.18) 

where A is the length of the killing vector. Also regular are a a = N a /X and g ab = ^(g a b + 
a a a b ) with 

^e# 3) -i, (4.19) 

9ab - S ab E H 3j5 . (4.20) 

It can be shown that g ab is a regular riemannian metric, at least in a neighbourhood of 
infinity. 

The next step is to realise that u a 

0Ji = -\ 2 e ljk D J a k (4.21) 

is curlfree, and can be replaced by a scalar u>, u a = D a u. Finally, Beig and Simon replace 
{\u)hj(j) A (A = 1,2) 

</>i = ^A- 1 (A 2 + ^ 2 -l), «/> 2 = 1a-V (4.22) 

The field equations now read 

V 2 <pA = 2R(f)A (4.23a) 
R ab = 2D a (j) A b h (j) A - D a ZD b Z, (4.236) 

with E = |(l + 40 A A )i 

Now we can iterate on (4.23) without difficulty and show that we can improve both 
the decay rates and the differentiability of (g ab , 4>a) essentially without limit. Eventually 
we reach high enough differentiability so that we can replace the weak derivatives with 
ordinary derivatives; we then can use Beig and Simon to prove analyticity. 

Thus we have proven: 
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Theorem. Given a solution to the Einstein equations which is asymptotically Hat , vacuum 
outside a region of compact support and possesses a timelike killing vector in the exterior 
region, such that g^ u —n^ u £ H^^iR 4 — B) for any 5 < in the natural coordinates (where 
d/dt is the killing vector) the metric must be analytic and Schwarzschildian at infinity 
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